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ABSTRACT 

The dense, clustered environment in which massive stars form can lead to interactions with neigh- 
boring stars. It has been hypothesized that collisions and mergers may contribute to the growth of 
the most massive stars. In this paper we extend the study of star-disk interactions to explore en- 
counters between a massive protostar and a less massive cluster sibling using the publicly available 
SPH code GADGET-2. Collisions do not occur in the parameter space studied, but the end state of 
many encounters is an eccentric binary with a semi-major axis ~ 100 AU. Disk material is sometimes 
captured by the impactor. Most encounters result in disruption and destruction of the initial disk, 
and periodic torquing of the remnant disk. We consider the effect of the changing orientation of the 
disk on an accretion driven jet, and the evolution of the systems in the presence of on-going accretion 
from the parent core. 

Subject headings: circumstellar matter — stars: formation 
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1. INTRODUCTION 

In the standard picture of star formation, stars form 
from the collapse and fragmentation of turbulent molec- 
ular cloud cores on a timescale of ~ 10 5 yr for individual 
stars to ~ 10 7 yr for large OB associations. Accretion 
onto the prestellar core from the surrounding envelope 
can occur spherically, or if the material has some angu- 
lar momentum, from the circumstellar disk. While the 
details of these processes are becoming understood for 
isolated low mass stars (M* < 10 Mq), the formation 
of clusters, multiple systems, and massive stars are less 
clear. Simply extending the process for low mass stars 
to massive star formation presents some problems. 

The Jeans mass of cloud cores containing massive stars 
is ~ 0.3 Mid, while obse rv ed stellar masses are ~ 50 
Mq HBonnell et alJll99l . IWolfire fc Cassinellil (fl987h 
argue that radiation pressure from the central star may 
halt or reverse spherical ac cretion for mas s es gre ater than 
about 20 M©. However. iMcKee fc T an (20030 suggest 
that high external pressure can overcome this obstacle. 
Calculations of accretion from th e envelope onto a disk, 
and from the disk onto the star ijYorke fc Bodenheimerl 
Il999t lYorke fc Sonnha lterll200 2[) have result ed in proto- 
stars of ~ 30M^. iKnrmholz et all (|2005albf l argue that 
anisotropics in the stellar radiation pressure and instabil- 
ities in the accretion flow can circumvent the radiation 
pressure problem. However, these scenarios deal solely 
with the formation of a star in isolation. 

Youn g massive stars a re usually found in dense clus- 
ters llCla rke et aLll2000D, often near the central region 
( Hillenbra nd fc H artman nlll998fl. and have a high multi- 
plicit y fraction ijStahler et all200GT:lGarcia fc Mermilliodl 
1200 1|) . These factors suggest that interactions between 
neigh boring sta r s are likely during cluster formation. 
Bon nell fc Bate) {2002) performed numerical simulations 
of accretion and merging in a cluster of 1000 identical 
mass stars. They found that 19 mergers occurred, driven 
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by the dynamics of small subgroups within the cluster. 
Though the stellar densities they required (~ 10 s star 
pc~ 3 ) are far in excess of those observed, it can be ar- 
gued that such high densities exist briefly during the 
initial collaps e of a high density molecu lar cloud core. 
iBonnell et all lfl99c| and IBonnell fc Bate! {2005) suggest 
a scenario wherein accretion onto a wide, low mass bi- 
nary leads to a tight binary with a system mass of 30-50 
Mq, which, with its larger cross section, could be driven 
to collision by an encounter with a third body. 

If one considers a cluster with number density 10 6 
pc~ 3 , a velocity dispersion of 3 km s _1 , and populated 
by typical pre-main-sequence protostars with masses 0.3 
Mq and radii 0.05 AU, one can estimate the mean time 
between collisions (Bin nev fc Tremainel 11987 !) as t co i ~ 
4 x 10 8 yr. This corresponds to a virialized cluster core 
with, e.g., 100 M (~ 330 protostars) and radius ~ 0.04 
pc. Even at such an extreme number density, this rate is 
not high enough for collisions to be a likel y star forma- 
tion mechanism. iBallv fc Zinneckerl l)2005|) suggest that 
circumstellar disks and envelopes could play a role be- 
yond accretion by increasing the interaction radius, and 
therefore assisting stellar capture. 

The cross section of a star-disk system is much larger 
than the cross section of a single star, with disks com- 
monly a fewxlO 2 AU in radius. Material ejected dur- 
ing an encounter carries away kinetic energy from the 
impactor, which expands the range of initial conditions 
that could lead to a merger. Additionally, by sweeping 
through the disk the impactor is slowed both by accret- 
ing material, and by the gravitational effects of wakes and 
spiral arms induced in the disk. Can these effects increase 
the merger radius of the protostars to ~ 100 AU? The 
increased gravitational focusing between a 20 Mq star in 
a mass-segregated cluster wit h density n+ = 10 5 p c~ 3 , 
similar to the density found in IBonnell et alJ l|2003li . can 
drive the timescale for encounters within 200 AU to less 
than 10 5 yr; encounters with stars more massive than 
1 Mq take place on a timescale ~ 10 5 yr (Moeckel & 
Bally 2006, in preparation). Can the ensuing interac- 
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tions lead to mergers, or to binaries which could evolve 
further through third body interactions or accretion onto 
the system? 

Several authors ha ve cons i dered e ncounters between 
star-disk sy stems. iHellerl (119951) . lLarwooc 
Bo ffin et all (|1998H and iPfalzner et all ll20oT 
body and SPH codes, while lHall et all 1)1996|) and 
lHalll l)1997|) used a reduced-three body approach to 
study encounters rangi ng from penetrating to distant. 
iClarke fc Pringlel l)1991[) approached the question of bi- 
nary formation due to star-disk encounters analytically, 
while lOstrikerl l)1994l) an alytically treat e d the ca se of 
more distant encounters. iWatkins et alJ (|1998albf ) per- 
formed SPH simulations of encounters wherein both stars 
have circumstellar disks. Of the numerical studies deal- 
ing with close encounte rs, none con sider stellar mass ra- 
tios m 2 /m 1 < 0.5, and lHall l(T997|) and iPfalznerl pOOl 
are the only studies that consider repeated encounters. 
In this paper we present the results of numerical simu- 
lations concerning repeated, penetrating encounters be- 
tween a massive protostar with a disk, and a less massive 
stellar impactor, to extend the parameter space consid- 
ered in previous studies toward a regime suitable for the 
early life of a high mass star in a star forming cloud core. 
The primary system is simulated in a late stage of forma- 
tion, when the ratio of disk to stellar mass M^/M* =0.1, 
though we consider one case with an accreting envelope 
surrounding the system. 

2. SIMULATION METHOD 

The large dynamic and spatial ranges inherent to 
the problem, combined with the lack of a conve- 
nient fixed geometry, point naturally toward the use 
of a Lagrangian method. We use the publicly avail- 
able smoothe d particle hy drodynamics (SPH) code 
GADGET-2 l|Snringe]l 1200.51) . The equations of mo- 
tion integrated by GADGE T-2 self-consistently include 
the so-called V ft terms ijSpringel Hz Hernauistl l2002t 
IMonaghanl I2002f) . Rather than advancing the internal 
energy, an entropic function is integrated, from which 
the thermal energy can be evaluated. Gravity is com- 
puted by means of a tree. 

GADGET-2 can accommodate a mix of collisional 
and collisionless particles. We model the stars as point 
masses interacting only through gravity, while the disk 
particles are subject to both gravitational and gas dy- 
namical forces. We modified GADG ET-2 so that st ellar 
particles behave as sink particles l|Bate et al.lll995() . ac- 
creting gas that falls within a set accretion radius. A 
particle is accreted only if it is energetically bound to 
the star, and if it has less angular momentum around the 
accretor than that of a Keplerian orbit at the accretion 
radius. Following accretion, the gas particle is removed 
from the computation. Its mass and linear momentum 
arc transfcrcd to the star, which is placed at the center 
of mass of the two particles. 

The accretion radius must satisfy two criteria. First, 
it should be small compared to the gravitational radius, 
r g = y/2GM acc /v 2 , with M acc the mass of t he acc retor 
and v its velocity through the gas. iRuffertl (|1096f) con- 
sidered the effect of an accretor's physical size relative 
to its gravitational radius in Bondi-Hoyle accretion. The 
physical radius of an accreting star is much less than it's 
gravitational radius; if r acc > 0.2r ff , the linear momen- 



tum transfer from directly accreting upstream material 
(versus from a gravitationally focused accretion stream) 
chang es sign. The accreting spheres considered in lRuffertl 
( 1996) are totally absorbing; our energy and angular mo- 
mentum tests help to some extent to decrease spurious 
accretion from upstream material. Ideally, r acc could 
be set strictly according to the criterion r acc /r g < 0.2. 
However, a resolution limit also needs to be taken into ac- 
count. If the SPH smoothing lengths of the gas particles 
nearest the accretor are large enough that they are neigh- 
bors with particles on the other side of the star, then they 
can be artificially supported from accretion by unphysi- 
cal gas pressure. Likewise, if the gravitational smoothing 
length, which we set to be approximately equal to the 
average SPH smoothing length, is larger than r acc then 
accretion will be stunted. 

A 5 M Q accreting object passing through a Keplerian 
disk around a 20 M star with a closest approach ~ 50 
AU will have r g > 25 AU for a prograde encounter, and 
r g > 5 AU for a retrograde encounter. In practice, for the 
prograde case and for inclinations such that the passage 
of the impactor through the disk is out of plane, r acc = 5 
AU is sufficiently small to avoid affecting the results. In 
retrograde encounters, this is clearly too large as r acc > 
r g . An accretion radius less than 1 AU is suspect with the 
resolution we have. Tests indicate that accretion radii in 
the range 0.75 AU< r acc < 1.5 AU produce comparable 
results; we choose r acc — 1.5 AU for simulations with 
inclination angles less than 45 degrees from retrograde. 

Since we are interested in encounters spanning multiple 
binary orbits, we require long term stability of the disks. 
Because of the artificial viscosity employed in SPH, disks 
have a te ndency to rapidly l ose mass onto sink particles. 
Following iBate et alJ l|1995D . we apply a corrective out- 
ward pressure force to particles close to the accretion ra- 
dius to compensate for the lack of pressure support from 
particles interior to the disk inner edge. This corrective 
force is only applied to particles in the disk around the 
primary star. The artificial viscosity formulation used in 
GADGET-2 is a variation of the standard M o naghan- 
Balsa ra viscosity ( Monaghan fc Gingoldl Il98fl iBalsaral 
1995), suggested bv IMonaghanl l)1997j) . utilizing a sig- 
nal velocity between two particles. This is implemented 
in GADGET-2 by 



siq 
V = Ci 



3wij, 



(1) 



where c is the particle's sound speed. The viscous tensor 
is then 



n y = 



siq 

a % w 
2 pij 



(2) 



Here w = Vij ■ tij/\T%j\ for approaching particles, and 
w = for diverging particles. Subscript ij denotes the 
arithmetic mean of the quantity for particles i and j, and 
a is the dimensionless viscous parameter. 

We have a dded one further mod i ficatio n to the viscos- 
ity, following IMorris fc Monaghan! (^597). Each particle 
has an associated a, evolving between the minimum a_ 
and maximum a+ via 



da 
~dt 



(3) 



Viscosity is thus kept at a_ unless the source term S 
acts upon it, after which it decays with a characteristic 
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Fig. 1. — Initial disk conditions. Left: Azimuthally averaged temperature, surface density, and Toomre Q profiles. Right: Face-on and 
edge-on log-surface density plots of the disk. Axes are in AU, and surface density is in Mq AU — 2 . 



time r = h/(kc), where k is a dimensionless parameter 
and h is the SPH smoothing length. The source term, 
meant to detect the presence of shocks, takes the form 

S = fa+max(-\\7 -v\,0), (4) 

with a+ our maximum allowed value of the viscous pa - 
rameter and / a shear viscosity limiter (iBalsaral 1995). 
We use the parameters a_ = 0.1, a + = 0.8 (typ- 
ical values for a in the unmodified viscosity imple- 
mented in GADGET-2 are - 0.75), and k = 2.5. Note 
that GADGET-2 defines the smoothing region over one 
smoothing length, so that this is equivalent to k = 5 for 
most other SPH codes. A simila r modification to t he 
viscosity in GADGET-2 is used bv lDolag et"afl <|2005|) . 

3. INITIAL CONDITIONS 

In order to study interactions between a star-disk sys- 
tem and an impactor, we first need a stable disk around 
the primary. We set up the disk with ~ 4.2 x 10 4 SPH 
particles, using the surface density profile 

£(r) = • (5) 

The vertical density structure at a given radius is given 

by 

p{r, z) = P°( r ) ex P[- 2 H(r) 2 J ' ^ 



where H(r) ~ r[c s (r) /v<j,(r)] is a temperature dependent 
scale-height and po the density at the mid-plane. The 
number of particles used is sufficient that the average 
smoothing length satisfies h < H(r)/2 for all but the 
innermost particles. The particles are initially placed in 
pairs, symmetric about the central star, at positions ran- 
domized slightly around a grid which would give the dis- 
tribution of equations J5J| and © . The simulation is run 
with just the primary star and the disk until the particle 
distribution reaches an equilibrium configuration. 
We use a temperature profile 

( r \- 1/2 

T{r) = T y—J , (7) 

where T is set depending on the mass of the star the 
particle is closest to. We implement a simple cooling 
scheme to maintain this temperature profile while the 
disk is unperturbed, 

du U- Ubase ,„n, 
~77 • (°J 

at t c 

The timescale t c is inversely proportional to the Keple- 
rian orbital frequency at the particle's radius, and Ub aS e 
is the internal energy corresponding to the temperature 
profile given by equation 10. The cooling time we use 
is short enough that the disk wou ld eventually fragment 
ljGammidl200S iRice et al.ll2005l) in the absence of the 
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Fig. 2. — Semi-major axis and eccentricity evolution for the various inclination angles. The top symbols are after the first encounter, 
followed by the values after subsequent even numbered encounters. The dotted line connects the values after 1 encounter, and the dashed 
line connects the values after 10 encounters. 



base temperature profile, below which the particles are 
not allowed to cool. Since we are interested in disks that 
are stable for long time periods, we choose initial condi- 
tions such that the Toomre parameter Q is greater than 
unity at all radii. Figure^ shows the temperature, den- 
sity and Toomre Q profiles of the equilibrium disk, as 
well as column density plots. 

4. ANALYTICALLY ADVANCING SIMULATIONS 

In order to save computational time, if an encounter 
results in a binary with a semi-major axis much larger 
than the disk radius the system is not integrated through 
apastron and back. Instead, once a simulation has ad- 
vanced so that the stars are far enough apart that tidal 
effects on the remnant disk are no longer important, and 
the binary's semi-major axis and eccentricity stop evolv- 
ing, numerical integration is stopped. This separation 
is typically several initial disk radii. By analytically ad- 
vancing the system along a Keplerian orbit until the ef- 
fects of the disk again become important, we can reduce 
the computation time by up to an order of magnitude for 
several of the cases considered. 

When a system is ready to be analytically advanced, 
the particles are grouped into one of three populations; 
those which have been ejected, those that are bound to 
the impactor, and those bound to the primary. If a par- 
ticle's total energy relative to the stars' center of mass is 
positive, it is ejected and removed from the calculation. 
The remaining particles are assigned to a star based on 
which potential well is deeper at that location. Though 
simplistic, this scheme gives good results once the stars 
separate enough that two clear populations are visible. 
At this point we treat the populations as point masses at 
their centers of mass, and determine an elliptical orbit. 
The two stars and their associated gas particles are ad- 
vanced analytically along this orbit until the point when 
the separation is 350 AU. 

Upon restarting, the particles in each group are placed 
using a circularizing scheme similar to that employed by 



iHalll l)1997j) . as follows. The angular momentum of the 
particle about the star is determined, as well as the sep- 
aration vector from the star to the particle. The sepa- 
ration vector is stretched to the distance of a Keplerian 
orbit with the same angular momentum. The velocity is 
made Keplerian, with the same orientation of the angu- 
lar momentum vector, and the particle is advanced along 
its new orbit for the length of time between stopping the 
simulations and restarting it. 

This circularizing scheme is an approximation, and one 
that is probably not accurate for material which is on the 
verge of being unbound, and is now on a highly eccen- 
tric orbit. However, most of the material is in a per- 
turbed disk close to the primary star, and circularizes 
fairly quickly if the simulation is allowed to continue. 
For this material, the procedure is reasonable as long as 
the orbital period of the binary is long compared to the 
orbital period of the disk particles. The viscous time 
scale for r > 20 AU is greater than the longest binary 
orbital period, so we can ignore viscous spreading of the 
disk between encounters. This technique is only applied 
to simulations with times between encounters ~ 10 3 yr; 
when the system evolves to the point where the time be- 
tween periastra is ~ 10 2 yr, comparable to the orbital 
timescale at the edge of the disk, numerical integration 
is allowed to continue uninterrupted. 

5. RESULTS 

In this study, we keep the mass of the primary star 
fixed at m p = 20 M Q , and the disk mass at md = 2 M Q . 
The main group of simulations are between that star- 
disk system and a TOj = 5 M Q impactor. The impactor 
is introduced on a slightly hyperbolic (uqo = 0.5 km s _1 ) 
orbit that, unaltered by the disk, would have a perias- 
tron distance of 50 AU. The initial separation between 
stars is 350 AU. This group of simulations covers the 
range of inclination angles i =[0,180] degrees, in steps of 
30 degrees, and each set of initial conditions is advanced 
through at least 10 encounters. The retrograde, coplanar 
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case (i — 180) was also run with an impactor mass of 1 
M© and a 50 AU periastron, and with a 5 M© impactor 
at a periastron of 100 AU. For simulations which are an- 
alytically advanced, the orbital characteristics discussed 
below are measured when the simulation is stopped; for 
those allowed to run continuously, they are measured at 
the binary apastron. 



5.1. Inclination Angle 

We introduce the 5 M© impactor such that the angu- 
lar momentum vector of the disk lies in the plane de- 
fined by the angular momentum of the orbit and the 
semi-major axis of the orbit, so that out of plane en- 
counters pass through the disk twice. An inclination 
of i = 0, a prograde encounter, is when the disk and 
orbital angular momenta are aligned; a retrograde en- 
counter, i = 180, has the angular momenta anti-aligned. 
All of the simulations considered here ar e pene t rating 
enc ounters, similar to SPH simulations bv lHellerl l)1995fl 
and Boffin et al]ll9 981. and reduced three-body calcula- 
tions bv lHall et alJ l)1996|) . though at a much higher mass 
range, and different mass ratios. Like those authors, we 
find that the dissipation of orbital energy is most efficient 
for retrograde encounters. 

Figurc|21shows the long term (over many binary orbits) 
evolution of the binary semi-major axis and eccentricity 
for each inclination considered. The dotted line connects 
the values after one encounter, and the dashed line con- 
nects the values after 10 binary orbits. The semi-major 
axis a is related to the specific orbital energy E and the 
total binary mass M = m p + m, by a = GM/2\E\. For 
the purposes of determining orbital characteristics, the 
gas mass is distributed between the primary and the im- 
pactor as described in section 0J Plotted in figure is a 
rough measure of the importance of disk disruption and 
heating by shocks in the dissipation of orbital energy. 
Disk disruption is measured by A(m p + nid), while the 
energy lost due to shock heating of the disk gas is mea- 
sured by the total cooling energy (see eq. [8]) calculated 
during the encounter. Disk disruption is maximal for a 
prograde encounter, while heating losses are maximal for 
retrograde encounters. 

Dis k disrupti o n is n ot a reliable tracer of energy trans- 
fer; lHall et all l|199fifl note that while a prograde en- 
counter is the most spatially destructive, a retrograde 
encounter has a larger energy transfer. Nonetheless it 
provides some idea about the mechanism behind the en- 
ergy loss from the orbit. An encounter with i — 150 has 
~ 1.4 times the energy loss from heating as a prograde 
or i = 30 encounter, but a larger semi-major axis after 
the collision. The dominant energy transfer mechanism 
for these nearly-prograde encounters is evidently not the 
shock heating of disk material but the gravitational in- 
teractions with, and disruption of, the initial disk, ev- 
idenced by over twice as much disk material being lost 
from the primary for i ^ 30 than for i = 15 0. This differ- 
ence in dissipation mech anism is noted bv iHellerl l)1995t) 
and iBoffin et alJ 1^998), and the general shape of the 
semi-major axis plot in figure [21 is consistent with those 
author's findings for \E or ut with inclination angle. 

r P eri = a(l - e), (9) 
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Fig. 3. — Disk disruption (solid) and energy loss from cooling 
(dash), relative to the maximum value for all inclinations. Disk 
disruption is a more important energy transfer mechanism for low 
inclinations, while losses from shock heating (measured by the sub- 
sequent cooling) dominate for retrograde encounters. 



with the eccentricity given by 

I 2E.P , . 

e =y 1+ 7PMi- (10) 

We run all simulations until cither the semi-major axis is 
essentially unchanged with further orbits, or the disk is 
destroyed, with two exceptions. The 90 and 120 degree 
inclination cases continue to evolve after 14 and 11 orbits 
respectively, and their disks retain ~ 1 M©. However, 
the elapsed times for these simulations are approaching 
3 x 10 4 yr. Note that the retrograde case has undergone 
27 orbits in less than 4 x 10 3 yr. For the environments 
we are considering, dense cluster cores, the timescale and 
size of the orbits (a ~ 1500 AU) for these cases are too 
large to ignore perturbations from other cluster members 
or the addition of mass from the cloud core. Plotted in 
figure^Jare the semi-major axes as functions of encounter 
number and time, and the half-mass radii and masses of 
the disks after each encounter. The disk's half-mass ra- 
dius is defined as the radius inside of which 50% of the 
gas energetically assigned to the primary is located. The 
general trend for all inclinations is a period lasting sev- 
eral encounters during which a steady decrease in orbital 
energy takes place. After a few encounters, this rate de- 
creases and the change in orbital parameters levels off. 
Note that the prograde case, while moderately effective 
at dissipating orbital energy over a few encounters, has 
the largest semi-major axis after 10 encounters. The rea- 
son for this is a combination of two effects. The high rate 
of disk disruption for low inclinations is a major factor; 
the prograde case has reduced the disk mass from 2.0 to 
0.5 M© after only 5 orbits, whereas the retrogra de case 
takes 10 orbits. This work and previous work (iHeilej 
[199 5 IHall et alJll99fit IBoffin et alJIT998t IPfalzner"etaL 
2005) verify the intuitive idea that with a larger disk 
mass, the orbital energy decay is greater. 

The evolution of the periastron distance also helps en- 
counters with a retrograde component achieve faster bi- 



6 





Fig. 4. — Left: Semi-major axis evolution for the various inclination angles, as a function of (top) encounter number and (bottom) time 
since the first encounter. For clarity, values in time are plotted as step functions measured at apastron. Top right: Disk half-mass radius as 
a function of encounter number. Bottom right: Remaining disk mass as a function of encounter number. Values are measured at apastron, 
or when the simulation is stopped for analytic advancing. 



nary decay during the disk destruction phase. Plotted 
in figure [S] are the periastra after each encounter for all 
inclinations. The periastra can be calculated from the 
orbital parameters by 

Thus, the periastron distance is a function of the or- 
bital energy and angular momentum. Note that E is 
negative for an elliptical orbit; for a given decrease in E, 
there is a critical change in J, above which the periastron 
will increase. Encounters with i < 90 don't lose enough 
orbital angular momentum (through accretion and gravi- 
tational interactions) to decrease the periastron distance, 
so while the disk radius is decreasing, the encounter dis- 
tance increases. This behavior is due in part to the lo- 
cation of the corotation r esonance region for the differ- 
ent prograde inclinations ijHall et alJll996|) . Material in 
the resonance region tends to t ransfer an gular momen- 
tum from the disk to the binary l)Ostrikerll994l) . working 



against the dissipative effects of accretion and material 
outside the corotation radius. 

Previous work on similar enco unters iHellerl Il995t 
IHall et alJlToM IBoffin et aLlflUM IPfalzner et all 12003) 
shows that with increasing periastron distance, the mag- 
nitude of the change in orbital energy decreases. Ret- 
rograde encounters, however, show a large decrease in 
periastra, down to approximately the disk half mass ra- 
dius. Comparing the disk characteristics for i = and 
i = 180 after 5 encounters (FigureEJ) with their periastra 
evolution, we see that while r$o ~ 22 AU for both cases, 
the closest approach for i = 180 is ~ 26 AU, still close 
to the disk, while for i = the closest approach is nearly 
60 AU, so that r per i > 2r so- 
la several cases with a prograde encounter, mass is 
transfered from the disk, resulting in a secondary disk 
around the impactor. This is illustrated in figure for 
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Fig. 5. — Periastron separation for each inclination as a function 
of encounter number. 



i = 30, the most striking example. In this case the 
disk is nearly orthogonal to the primary disk, a sce- 
nario tantalizingly similar to the non-coplanar disks seen 
in young multiple sy stems such as HH111 and HH121 
ijReipurth et alJlT999j) . The disk survives (or is replen- 
ished) during subsequent encounters. Increasing perias- 
tron separation and rapid disk destruction are seen in 
both the prograde encounter and for i = 30; we suspect 
that the increased interaction between the disks causes 
the out-of-plane encounter to end up in a tighter binary 
An effect noted in pre vious work llHerlerll995tlLarwoodl 
119971 iBoffin et alj fl998) is that out-of-plane encounters 
lead to tilting of the remnant disk. After several encoun- 
ters, this effect is seen even in the in-plane, retrograde 
case. When followed through several orbits, this effect 
manifests itself as stepped precession, shown in figure [7| 
The orientation of the disk is measured as the average 
angular momentum of gas particles about the primary 
star. In figure we consider only the inner disk, r < 30 
AU. The disk undergoes a slow precession during most of 
the orbit, with periodic jumps associated with periastron 
passages. There is warping present in the disks, shown in 
figure |H1 for the 30 degree inclination case. Plotted is the 
difference in the orientation of annuli 10 AU wide at four 
different times, relative to the central annulus. This an- 
gle is measured by i — cos \L a -Lb/dLaHLbl)], with L a 
and Lfc the summed angular momenta of all particles in 
two annuli. If the orientation of an outflow is tied to the 
inner disk angular momentum vector , this is potentially 
an observable effect l|Bate et al.l2 000). Shown in figure^ 
is a crude estimation of what such an outflow might look 
like, again for the retrograde case. We assume that a 
jet is launched along the disk angular momentum vector, 
at a uniform velocity for the entirety of the simulation. 
The units are arbitrary and depend on the velocity of 
the outflow; for a 10 km s^ 1 jet, the long axis of the jet 
extends ~ 8000 AU. It isn't until ~ 2000 years, or 10 
binary orbits, that the jet begins to precess noticeably. 
As the encounters become more frequent and the disk 
becomes slightly non-coplanar to the binary orbit, the 



precession becomes more pronounced. The sharp kinks 
in the jet are a result of the torque being applied almost 
totally during periastra passages. 

5.2. Evolution of a Binary in an Accreting Envelope 

Material from the parent cloud that accretes onto a 
binary can alter the binary orbital elements and mass 
ratio, effects that have been studied by several authors 
(lArtymowicd Il983t iBatej Il997t iBate k Bonnelil Il997t 
iBatel 120001) . lArtvmowica l)1983t) suggests that low an- 
gular momentum gas falling nearly radially onto the sys- 
tem will result in a decrease of the semi-major axis, a 
result confirmed bv IBate k Bonnelil l)1997j) . IBatel (|2000f) 
shows that the decrease in separation can reach an order 
of magnitude or more, provided that the mass accreted 
onto t he system i s sever al to 10 times the total stellar 
mass. IBate et all 1)2002(1 describe a system whose com- 
ponents increase their masses several fol d, while the sepa- 
ration decreases from ~ 100 to ~ 10 AU. lBonnell k Bate! 
(2005) show that accretion onto a low mass (m ~ 1 M Q ) 
binary with an initial separation ~ 1000 AU can lead 
to a close (separations ~ 1 AU or less) binary with a 
total mass ~ 30 — 50 M a . In all of these studies the 
accreted mass is comparable to, or orders of magnitude 
larger than, the initial system mass. 

In order to see a comparable effect with a total stellar 
mass of ~ 25 M Q , a much larger mass must be added 
to the system. If the system mass is to double in less 
than 10 5 yr, an average accretion rate of M > 2.5 x 10 -4 
Mq yr -1 is required, which is fairly high, though perhaps 
not unreasonable for the formation site of a 20 M star 
ijMcKee k Tarjl2003l) . Given the previous results, even 
this high rate would result in a hard binary only after 
~ 10 5 yr, which is long compared to the encounter time. 
It seems that in massive systems such as those considered 
here, accretion from a surrounding envelope will not be a 
large factor in hardening binaries. However, the previous 
work has been mostly concerned with less massive, and 
less eccentric binaries. In order to explore this issue, we 
consider the effect of a massive accreting envelope in the 
retrograde case. 

We add an in-falling envelope with M = 10~ 3 
M Q yr -1 , consisting of ~ 5 x 10 4 particles, originating 
at a radius of 500 AU from the center of mass of the sys- 
tem. The in-fall is spherically symmetric, purely radial, 
and is introduced with a velocity equal to the free fall 
velocity at 500 AU. The lack of rotation is unrealistic, 
but chosen for maximal effect on the orbit. The cutoff 
radius of 500 AU was chosen to keep the simulation time 
reasonable. The effect of any ejected mass on shaping 
the in-fall at larger radii is thus effectively negated. A 



in-fall morphology significant!) 


r lYorke k Bodenheimerl 


1999tlYorke k Sonnhalterll2002 


iKrumholz et all2005b|) . 



Spherical in-fall, and thus the lack of angular momen- 
tum in the envelope, should lead to a maximal decrease 
in semi-major axis, so this simulation is illustrative of 
the most extreme case. 

Plotted in figure ITUl are the semi-major axis and mass 
evolution of the retrograde case with and without an in- 
falling envelope. Spikes in the semi-major axis plot are 
an artifact of our orbital parameter calculation scheme 
during periastra approaches. The gas mass plotted is 
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Edge-on snapshots of the first encounter for i = 30. Plotted is log column density, in units of M© AU 
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Fig. 7. — Top: Inner disk orientation (r < 30 AU) as a function 
of time, for i = 180. Plotted are the azimuthal angle 9 (solid) and 
ten times the polar angle (dashed) relative to the initial disk 
orientation. Bottom: Binary separation for the same simulation. 
Pcriastra are marked by dotted lines. 



that which is energetically bound to the primary star. 
The initial binary decay is faster for the case with the 
in-falling envelope, though the final state of the system 
is similar to the case with no envelope. Comparing the 
masses of the stars for the two cases, it is apparent that 
most of the infalling mass is accreting onto the central 
star. At t ~ 3000 yr the impactor in the envelope case 
has a mass rrii ~ 5.4 Mq, only slightly more than rrii ~ 
5.3 Mq in the absence of an envelope. 

At t ~ 3000 the primary mass without an envelope 
has leveled off at m p ~ 21.1 Mq, and the disk mass (the 



difference between 



and 



L gas 



) has decreased to 



ma ~ 0.1 Mq. With the addition of a spherical infall, 
the disk mass is minimally effected, again approaching 
negligible mass after t ~ 2500 yr, at which point the pri- 
mary star enters a period of steady growth at M ~ 10~ 3 
Mq yr -1 . The primary itself has grown to m p > 23 M Q ; 
it appears that the infall is being tunneled almost en- 
tirely onto the primary, spending little time in the disk. 
Because we were searching for a purely accretional ef- 
fect as opposed to the influence of a circumbinary disk 
l)Bate fc Bonnelllll997j) . the infall was set up with no 
angular momentum. In this case, once the system has 
evolved to a near steady state, with grazing disk encoun- 
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At = 218 yr 
• t. = 265 yr 
* t = 313 yr 




Fig. 8. — Warping of the 30 degree inclination disk, measured as 
the difference in orientation between annuli 10 AU wide and the 
central annulus. Times are measured from the first periastron. 
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Fig. 9. — Top: A toy jet launched from the same disk as figure 
171 The jet is assumed to launch perpendicular to the inner disk, 
and with a uniform velocity throughout the evolution. 



ters slowly eroding the binary energy and circularizing 
the orbit, the stellar masses dominate the mass accreted 
from the envelope and the system is minimally affected, 
as expected. 

The total mass accreted from the envelope over 3000 
yr, while small compared to the stellar masses, is greater 
than that of the initial disk; if the angular momentum of 
the envelope is sufficient to accrete into a disk around the 
primary at a radius comparable to the binary orbit, the 
initial period of rapid orbital decay may be extended and 
tighter binaries could be formed. However, accretion of 
low angular momentum gas from the parent core seems 
unlikely to significantly alter the orbits of the massive 
capture formed binaries studied here. 

6. SUMMARY AND CONCLUSIONS 

We have presented a suite of simulations which con- 
sider the dissipative effects of a circumstellar disk on a 
capture formed binary, in the context of a massive pro- 
tostar encountering a less massive cluster member. The 



systems were followed for 10 or more orbits, in many 
cases leading to almost total destruction of the original 
disk. We find that for the parameters studied, the in- 
teractions result in a soft, eccentric binary system. The 
energy change associated with different inclination angles 
is in agreement with previous simulation results. The pe- 
riodic torquing of the disk by the passage of the impactor 
may be observable in the orientation of a jet or outflow. 

These simulations apply to a massive protostar that 
has accreted most of its final mass. The disk and enve- 
lope mass is ~ 10 % of the stellar mass, analogous to the 
class I phase of a low mass protostar's evolution. Previ- 
ous studies of stellar encounte rs with star-disk systems 
l)Hellerfll995tlBoffin et alll998j) concluded that disk cap- 
ture rates were too low to contribute significantly to the 
observed multiplicity of low mass stars. All but the least 
dissipative encounters simulated here result in a binary 
with a semi-major axis < 500 AU within 10 4 yr. Extend- 
ing the analysis of binary capture rates to the regime 
where the disk mass and impactor masses are compara- 
ble, and much less than the primary, may show that the 
high multiplicity of massive stars such as the Trapezium 
stars in Orion could result from disk encounters in the 
late stages of the cluster's formation. 

In order to follow the encounters through several or- 
bits, the parameter space covered here is small; only a 
single impactor mass and initial periastron were consid- 
ered. The mass of the impactor was chosen as the mini- 
mum such that a merger would lead to significant growth. 
As the impactor's mass approaches the primary's, disk 
destruction occurs over fewer orbits and less time is avail- 
able to harden the binary. The apastra of the resultant 
binaries are on the order of 100s of AU, and the periastra 
of the most dissipative cases remain above ~ 25 AU. For 
the parameters studied, it seems that in the absence of 
further effects, disk dissipation alone cannot lead directly 
to mergers. 

Accretion onto the system from the parent core is un- 
likely to lead to significant hardening of the binary over 
any short timescale, as the stellar masses considered are 
too high for reasonable mass accretion rates to have a 
large effect. Thus if disk mediated collisions are to con- 
tribute to the growth of a massive star, they must occur 
earlier in the cluster's evolution, when the disk and en- 
velope are more massive and perhaps more dissipative. 
Even if this is the case, there must be an efficient disk 
replenishment mechanism in order for more than a single 
merger to occur. However, at least two further mecha- 
nisms seem likely to come into play. 

The disks considered in this work are Toomre stable; if 
a colder disk is considered, gravitational fragmentation 
within the disk will change the nature of the interac- 
tions, potentially leading to a more rapid decrease in the 
binary separation as smaller bodies are ejected from the 
system. In a similar dynamical effect, a binary with a 
semi-major axis ~ 100 AU presents a large target for 
interactions with other encounter partners in the dense 
cluster core. In this scenario the capture and gradual 
tightening of the initial binary, modeled here, provides 
the initial conditions for encounters with other cluster 
stars. The dynamics of n > 2 bodies plus the remnants 
of a disk are beyond the scope of this work, but are an 
enticing direction for further study of systems such as 
these. 
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Fig. 10. — Left: Semi-major axis evolution with time for i=180, with (solid) and without (dot) an accreting envelope. Apastra are 
marked with symbols. Right: Impactor and primary masses, and mass of gas associated with the primary for the two cases. 



In summary, it seems that disk interactions in the late 
phase of protostellar evolution considered here do not 
provide a direct path to merger. Further accretion from 
the natal cloud will not be enough to induce collisions 
from the resultant binaries. In order for disks such as 
those modeled here to play a role in any collisional growth 
of massive stars, the interactions must involve third bod- 
ies acting on the capture formed binary, through either 
induced fragmentation of the disk or from other cluster 
members. Another, yet unmodeled possibility is that the 
encounters leading to a merger take place during an evo- 



lutionary phase more like a class protostar, with the 
mass of the primary comparable to that of the disk sur- 
rounding it. Encounters with such systems have only 
been considered with impactors of similar mass to the 
primary, and not in the limit where the disk mass is 
greater than the impactor. 

We thank the referee, Ian Bonnell, whose comments 
helped improved the manuscript. This work was sup- 
ported by the NASA Astrobiology Institute through a 
Cooperative Agreement with the University of Colorado, 
through grant NNA04CC11A. 
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